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ABSTRACT 

There is now strong evidence suggesting that the ^^CO J=l-0 transition, 
widely used to trace H2 gas, significantly underestimates its mass in metal-poor 
regions. In spiral disks such regions are found in large galactocentric distances 
where we show that any unaccounted H2 gas phase is likely to be diffuse (n ~ 
5-20 cm-3) and warmer (Tkm ~ 50 - 100 K) than the cool (Tkm ~ 15-20 K) 
CO-luminous one. Moreover we find that a high value of the H2 formation rate on 
grains, suggested by recent observational work, can compensate for the reduction 
of the available grain surface in the metal-poor part of typical galactic disks and 
thus enhance this CO-poor H2 component which may be contributing significantly 
to the mass and pressure of spiral disks beyond their optical radius. 

Subject headings: dark matter: galaxies — interstellar matter: molecular gas — 
galaxies: individual: NGC 891, Galaxy — molecules: CO, H2 

1. Introduction 

The use of ^^CO, the most abundant molecule after H2 itself, to trace H2 gas is a widely 
used practice since the discovery of bright ^^CO J=l-0 emission in Orion (Wilson, Jefferts, 
& Penzias 1970). The most easily excited H2 transition, the quadrupole S(0): Ju — Ji = 2 — 0, 
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corresponds to AE2o/k ~ 510 K, much too high to be significantly populated in the cool ISM 
(Tkin ~ 10-50 K) revealed by the ^^CO J=l-0 transition (AEio/k ~ 5.5 K). Furthermore the 
large optical depths of the latter (rio ~ 5 — 10) reduce its critical density of n^^^ ~ 380 cm~^ 
(for Tkin = 50 K; Kamp & Zadclhoff 2001) to - n^r°^(l -e-^io)r{o^ ~ 40 -80 cm-^ similar to 
the lowest average densities observed in Giant Molecular Clouds (GMCs). Hence the ^^CO 
J=l-0 transition is expected to be excited even in the coldest and most diffuse regions of 
the molecular ISM. 

However its large optical depths (e.g. Martin, Sanders, & Hills 1984; Sage & Isbell 
1991; Falgarone 1998), do not allow a straightforward use of its luminosity together with 
a [^^C0/H2] abundance ratio for H2 mass estimates. The widely used Xco = M(H2)/Lco 
factor (where Leo is the velocity/area- integrated brightness temperature of ^^CO J=l-0) is 
based on the notion of an effectively optically thin medium, where the ^^CO line emission 
arises from an ensemble of small, radiatively de-coupled cells not overlapping in space or 
velocity. Then the high optical depths arise locally within the individiial cells and the 
observed wide line profiles result from their macroscopic motions rather than their intrinsic, 
much narrower linewidths (e.g. Martin ct al. 1984; Young & Scoville 1991 and references 
therein). This picture, along with the observed virialization of molecular clouds (see e.g. 
Larson 1981) over a wide range of scales (~ 0.1 — 100 pc) allows a statistical derivation of 
Xco as an ensemble-average that is relatively insensitive to local molecular cloud conditions 
(e.g. Dickman, Snell, & Schloerb 1986; Young & Scoville 1991; Bryant & Scoville 1996) and, 
for scales > 10 pc, is expected to yield the correct H2 mass to within a factor of two. 

The existence of such cells is supported by high angular and velocity resolution studies 
of galactic clouds (e.g. Falgarone et al. 1998; Tauber 1996) while the statistical robustness 
of Xco for a variety of conditions has been verified using both observational (e.g. Young & 
Scoville 1991 and references therein), and theoretical methods (e.g. Kutner & Leung 1985; 
Dickman, Snell & Schloerb 1986; Maloney & Black 1988; Wolfire, HoUenbach, & Tielens 
1993; Sakamoto 1996; Bryant & Scoville 1996). However these studies also suggest a likely 
failure of this method in two particular environments, namely a) galactic and/or starburst 
nuclei where non-virial motions cause the luminous ^^CO emission to overestimate the H2 
mass (e.g. Dahmen et al. 1998; Downes & Solomon 1998), b) metal-poor regions where 
it seriously underestimates the H2 mass (e.g. Maloney & Black 1988; Arimoto, Sofue, & 
Tsujimoto 1996; Israel 1988, 1997, 1999). The last case is the most serious one since then 
^^CO emission is usually faint and thus precludes the observations of more optically thin 
isotopes hke ^^CO, C^^O that could yield the H2 mass without using the Xco factor. 

In this paper we show that the decreasing metallicity with galactocentric radius in disks 
and a larger H2 formation rate open up the possibility of large amounts of H2 gas being 
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there hitherto undetected by the conventional method. This phase will be warm, diffuse and 
devoid of CO, and may contribute significantly to the pressure and total mass of the disk. 

2. Spiral disks: H2 formation in the metal-poor parts 

Molecular hydrogen in interstellar conditions forms from HI association onto dust grains 
(e.g. Gould & Salpeter 1963) and deep sub-mm imaging shows the grain distribution extend- 
ing beyond the optical and well into the HI disk (Nelson, Zaritsky, & Cutri 1998; Xylouris 

et al. 1999), and possibly out to ~ 2R25 (Thomas et al. 2001, 2002; Guillandre et al. 2001). 
Thus the basic "ingredients" for H2 formation can be found to the outermost parts of the 
disk as defined by the HI distribution itself. Using the results of Federman, Glassgold, & 
Kwan (1979), with an exponent of m=1.5 for the H2 self-shielding function, we find that an 
HI cloud of radius Rd (pc), average volume density n and temperature T, illuminated by a 
FUV field strength Go (Go(solar)=l) starts turning molecular when s > 1, where 




The density Umin is the minimum required for the onset of an HI-to-H2 transition and z = 
Z/Zq is the normalized metallicity. We kept the explicit temperature dependence of the H2 
formation rate Rf = SfT^/^, assumed that it scales linearly with total grain surface and thus 
metallicity, and included the factor /x = Sf/Sf°\ to parametrize for a value different than the 
currently adopted Sj°^ = 3 x 10"^^ cm=^ s"^ K-V2 (Jura 1975a). 

Over an HI disk it is T(HI) ~ 80 — 150 K, with a slight increase with galactocentric 
distance Rgai (e.g. Braun 1997), a temperature range typical of the warm neutral halos 
around molecular clouds (e.g. Andersson & Wannier 1993) where the HI-H2 transition zone 
lies. Thus no significant variation of s is expected because of a strongly varying T(Rga,i). 
The infiuence of the z/Go ratio is less straightforward since significant metallicity gradients 
exist in spiral disks (e.g. Henry 1998; Garnett 1998 and references therein) along with the 
well known exponential optical (and hence FUV) light profiles. Their eff'ects have been 
thoroughly modeled by Elmegreen (1989, 1993) and applied to spiral galaxies by Honma, 
Sofue, & Arimoto (1995). The latter find the metallicity gradient to be the main cause of 
the so called molecular front, the steep H2-to-HI transition at a particular Rtr deduced for 
most spirals from CO imaging. Evidently a larger H2 formation rate will alter this picture 
by "pushing" this front to encompass a larger portion of a typical disk. Rates that are 
~ 4 — 8 times higher than the standard value are compatible with observationally deduced 
upper limits (Jura 1974), since they always contain an n — Sf "degeneracy" that can bias 
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Sf estimates towards lower values (e.g. Jura 1975a,b; van Dishoeck & Black 1986). Recent 
results from ISO observations of H2 rotational lines from photodissociation regions (Habart 
et al. 2000; Li et al. 2002) strongly favor n = Sf/Sf^ > 5 (see also Sternberg 1988 for early 
indications of such high values), which will raise the value of s at all Rgai. 

This rescaling has no effect for those regions of the disk where the old values are already 
s > 1 and which CO imaging shows them to be indeed H2-rich, but now a larger portion of 
the disk is expected to have s > 1 and thus conditions favorable for an Hl-to-H2 transition. 
A rough estimate of this new radius can be obtained from the results of Honma et al. (1995) 
that show this transition to occur very sharply at a particular Rtr- The expected FUV 
volume cmissivity jo, its dust absorption coefficient kd, and the metallicity follow a similar 
spatial distribution (Honma et al. 1995) hence, neglecting the H2 self-shielding contributions 
to the total absorption coefficient corresponds to z/Go = z/(jo/k(j) oc z. For an exponential 
disk profile and s(Rtr) ~ 1 it is 

Rtr(/i) = Rtr(l) + (ln/x)Re, (2) 

where Re is the scale-length of the exponential distribution of optical light and metallicity. 
Thus the molecular front will now reside at least (ln5)Rc ~ 1.6Rc further, past that inferred 
by CO imaging. In typical disks the later is at Rtr(l) = (0.5 — l)Rc (e.g. Young & Scoville 
1991; Regan et al. 2001), thus Rtr(Ai) ~ (2 — 2.5)Re places the molecular front well inside 
a typical HI disk. A flattening of the abundance gradient, and thus of dust grain surface 
per H nuclei, at large Rgai (Henry 1998) and H2 self-shielding will only strengthen the above 
arguments by yielding still larger values of s in the outer parts of the disk. 

Therefore it now seems that the molecular front in spirals, deduced from CO imaging, is 
in reality a (low) metallicity- induced ^^CO detection threshold, and the necessary conditions 
for H2 formation exist over a much larger part of a typical spiral disk. This is not surprising 
since it is known that Xco oc — z~^ (e.g. Arnault et al. 1988; Israel 1997), which alone 
implies a Rtr larger than that found by using ^^CO imaging and a constant value of Xco- 

2.1. The role of pressure in the HI/H2 transition 

In obtaining the previous result we assumed no significant changes of (n, R^) of the 
various clouds with Rgai, which leaves only the z/Go variable to define Rtr. In a forthcoming 
paper we will present a more detailed analysis of the HI-H2 transition without such sim- 
plifications, yet a re-formulation of s in terms of cloud mass M and boundary pressure Pe 
suggests that the aforementioned assumptions do not subtract from the generality of our 



- 5 - 



result. Using the results of Elmegreen (1989) (his S function and s are related simply as 
s — S^/^), we obtain 

^ ^^'^[w cm-^Kj \10^Mq) VIOOkJ VGoJ ' 

It is clear that unless M is a strong function of Rgai, a radially dechning Pe is the only other 
important factor besides /iz/Go in determining Rtr. The increase of s with decreasing M may 
seem counterintuitive but it is due to the fact that small clouds have higher average densities 
than big ones with the same boundary pressure. This stems from the widespread virialization 
of most of the clouds, which in combination with the observed velocity dispersion vs cloud 
size relation (e.g. Larson 1981), yields n oc R^" [a ~ 1). In the metal-rich part of a typical 
spiral, where CO can be relied upon as an II2 gas mass tracer, the bulk of the H2 mass is 
distributed in a cloud hierarchy with M ~ 10^ M0 being typically the upper limit. The 
value of the latter is not observed to change with galactocentric distance (e.g. Elmegreen 
& Elmegreen 1987) except in galactic centers where tidal effects may truncate such a cloud 
hierarchy at higher minimum densities and cloud masses (e.g. Maloney & Black 1988). 

The measurement of Pe, especially as a function of Rgai, is not easy and deducing 
it by assuming equipartition with the random cloud motions (e.g. Honma et al. 1995) 
relies on a known total gas surface density, which is what needs to be determined. This 
task is easier for the Galaxy whose disk we consider typical in that respect. In this case 
a pressure probe sampling its outer parts, namely diffuse CO-luminous molecular clouds, 
suggests pressures Pe(12kpc) ~ (1 — 2) x lO'' cm^^ K (Heyer, Carpenter, & Snell 2001), which 
are higher than expected and similar to the Solar neighborhood value. Further out pressures 
of Pc(17kpc) ~ 5 X 10'^cm~^ K are typical, with few clouds reaching up to ~ 10^ cm~^ K 
(Brand, & Wouterloot 1995). It is straightforward to show that since 0.25 ^ z/Go ^ 1 (at 
Rgai ^ 17 kpc), for II = 5, such Pe values yield s > 1. Thus the H2 formation conditions 
are fulfilled out to ~ 2.5Re(~ 17 kpc). This radius is well past the location of the ^^CO 
J=l-0 brightness peak at ~ 5 kpc or any significant ^^CO emission (e.g. Blitz 1997) and 
well into the HI disk. Interestingly HII regions, another potential pressure probe, may also 
hint at larger than expected pressures in the outer Galaxy (Rudolph et al. 1996) and dwarf 
irregulars (Elmegreen, & Hunter 2000). In both of these cases low metallicities make ^^GO 
a poor indicator of H2 gas which may thus be responsible for the high pressures. 
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3. Moleculsir gas at Icirge galactocentric radii: CO-deficient and wairm 

The arguments for favorable conditions for an HI-to-H2 phase transition presented here 
do not and could not yield a mass estimate for the associated H2 gas. Nevertheless the 
simple fact that this gas resides untraceable by CO in the metal-poor parts of the disk can 
offer significant insights regarding its state. An obvious consequence of the low metallicities 
is reduced dust shielding which allows even the weaker FUV radiation field at large Rgai's 
to dissociate the mainly dust-shielded CO and confine it to much smaller volumes while 
leaving the largely self-shielding H2 intact. A typical range of Go ~ 0.3 — 0.6 and an ~ 5 
yields a minimum H2 formation density of Uf ~ 5 — 7 cm~^ for cloud sizes of ~ 10 — 50 pc. 
Observations of CO-luminous molecular clouds yield average densities of ~ 50 — 500 cm^'^ 
irrespective of galactocentric distance, hence n ~ 5 — 500 cm~^ and z ~ 0.15 are appropriate 
inputs for Photodissociation Region models (PDR) of the molecular gas at large Rgai. Similar 
modeling of metal-poor gas has been performed before but at much higher UV fields and 
gas densities than here (e.g. Pak et al. 1998). 

We investigate the gas abundance of warm H2 using a PDR code that models the 
chemistry of a diffuse cloud represented as a semi-infinite slab, with a pseudo-time-dependent 
multi-point chemical simulation (Viti et al. 2001). The UMIST chemical database (Millar et 
al. 1997) was used with some changes (see Viti et al. 2001). The heating, the cooling and CO 
emission is as described by Taylor et al (1993); we have however made the following changes: 
i) the thermal balance is achieved using Ridder's method, a fast non-linear convergence 
algorithm, ii) the heating mechanisms have been revised and now include also large grains 
and PAHs (Weingartner & Draine 2001), as well as UV pumping (not important for low Go's 
and densities) and chemical heating, iii) the collisional rates have been updated to include 
all the dominant collisional partners for such an ISM phase (H, H2, He). All our runs were 
performed assuming a microturbulent line formation with a typical width of ~3 km s~^, and 
following Tielens & HoUenbach 1985 the level populations are updated from the outside in. 
The heating caused by the decay of the turbulent gas is negligible. 

The results are summarized in Figure 1 where it becomes apparent that in diffuse gas 
of sub-solar metallicity the CO-deficient H2 gas dominates the cloud. This radiation-driven 
spatial segregation also maintains a temperature difference in which the pure H2 phase stays 
always warmer than the CO-luminous one. This simply reflects the different dominant 
coolants present, namely ^^CO and its isotopes for the CO-luminous phase, and CII, H2 
for the CO-deficient one. Average temperatm^es of Tk ~ 15 — 20 K are observed for most 
quiescent CO-luminous molecular clouds in the disk (with T^ ~ 30 — 60 K found only in 
star-forming warm "spots") and remain almost constant with Rgai- In the Galaxy this is 
beautifully demonstrated by COBE observations of CO J-l-1— (J=0-7) lines that show 
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invariant relative strengths suggesting Tk ~ 20 K over the entire disk except the Galactic 
center (Fixsen, Bennett, & Mather 1999). This temperature is similar to those of several 
individual molecular clouds, some as far as Rgai ~ 30 kpc (Digel, de Geus, & Thaddeus 1994), 
and also similar to that measured by COBEioi the bulk of the dust in the Galaxy (Sodroski 
et al. 1994; Lagache et al. 1998) and by sub-mm imaging in other spirals (e.g. Alton ct al. 
1998). Such gas temperatures are indeed expected for the dust-shielded (Ay > 4) CO-rich 
H2 clouds, where CO cooling (Goldsmith, & Langer 1978; Neufeld, Lepp, & Melnick 1995) 
lowers Tk to that of the concomitant dust whose temperature is weakly dependent to the 
ambient ISRF and thus to Rgai. 

For the unshielded regions of the molecular ISM, now expected to contain the bulk of 
the H2 gas mass in the metal-poor parts of spiral disks, lack of the CO coolant will keep 
Tk > Tdust (e.g. Tielens, & HoUenbach 1985). Its temperature range can be found from 



where the coohng (A) due to the CII and H2 lines is balanced by the photoelectric heat- 
ing from dust grains Fph (e.g. Tielens & Hollenbach 1985) and cosmic rays. Unlike solar 
metallicity environments, the cooling due to H2 lines may now be important and is included. 

The CII emission emanates from a 2-level system, hence following Hollenbach, Taka- 
hashi, & Tielens (1991) and for optically thin emission it is 



where n^ ~ 3 x 10^ cm ^ is the transition's critical density, and all carbon is assumed in the 
form of CII, mixed with a fully molecular gas phase at an abundance of [C/H] = 3 x 10~^ z. 

For the expected range of densities and Tk ^ 200 K, (a typical upper limit for the Cold 
Neutral Medium (CNM) in spiral disks), only the S(0) and S(l) lines rotational lines of H2 
can be appreciably excited and thus contribute to the cooling. Moreover, since no radiative 
or coUisional processes "mix" ortho-H2 ( J=odd) and para-H2 ( J=even) states we can consider 
S(0) and S(l) practically arising from two different molecules, each approximated by a 2-level 
system, therefore 



Acn + Ah2 



(4) 





Ah, = 2.06 X 10 
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where lop = (n3 + ni) / (n2 + no) is the ortho-to-para H2 ratio, ncr(O) ~ 54 cm ^ is the critical 
density of the S(0) transition and Rio = As(i)/As(o), namely 



The critical density of the S(l) transition is ncr(l) ~ 10^ cm ^, and this term dominates for 
the warmest and densest gas in the chosen parameter space. 

We adopt the latest photoelectric heating model of Weingartner & Draine (2001) with 
Fph proportionally scaled by the normalized metallicity factor z and parameters correspond- 
ing to carbonaceous grains (the last two models in their Table 2). The assumed electron 
abundance is the average for the CNM, under similar set of (z,Go), i.e. ne/2n ~ 5 x 10~^ 
(Wolfire et al. 1995). Finally we include a heating factor from cosmic rays Fcr = 6.4 x lO^^^n 
ergs s~^ (Goldsmith & Langer 1978), which is < 10% of the photelectric term. 

For Go = 0.3 — 0.6, z = 0.15 and rop = 1 — 3, both a diffuse (~ 5 — 20 cm~^) and a dense 
(~ 100—500 cm~^) H2 gas phase satisfy (4) with the former being warmer (~ 50—100 K) than 
the latter (~ 20 — 30 K). At least for the Galaxy, these solutions can be further constrained 
by the values of the total pressure Pg < (5 — 20) x 10^ cm~^ K, at Rgai > 8 kpc, where 



(e.g. van Dishocck & Black 1986). The second term expresses the non-thermal part of 
the pressure due to the macroscopic motions of the gas, which are assumed of turbulent 
origin with being the one-dimensional FWHM of a gaussian distribution. The velocity 
dispersions of CO-bright molecular clouds or HI in face-on spirals yield cr^; ~ 5 km s~^ 
(Combes 1999 and references therein) hence the aforementioned constraints on the total 
disk pressure correspond to n ~ 5 — 15 cm~^ (for T < 100 K). Thus a metal-poor, warm and 
diffuse H2 gas in the outer parts of spiral disks is likely to exist, and may indeed be the type 
of gas detected through its S(0) line with ISO in the edge-on spiral NGC 891 (Valentijn & 
van der Wcrf 1999), with a reported temperature of T ~ 90 K and a mass that outweighs 
that of HI by factors of ~ 5 — 15. 

Several issues regarding such a gas phase remain open, namely its thermal stability, 
low star formation rates, its exact relation to the HI gas distribution and the dynamical 
evolution of such potentially massive gaseous disks. Moreover an enhanced H2 formation 
rate, if indeed at work, may alter relevant aspects of the "standard" ISM picture (e.g. gas 





Pe = nkeT + 56.8n(7,^ 



(8) 
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heating) that must be taken into account in a detailed analysis. These and other issues will 
be the subject of a future paper. 

4. Conclusions 

In the light of mounting evidence that the standard method of using rotational line 
emission of the ^^CO molecule to trace H2 underestimates its mass at low metallicities we 
examined the state of such gas in the disks of spiral galaxies at large galactocentric distances. 
Our results can be summarized as follows, 

A high H2 formation rate along with the presence of dust throughout a typical HI disk, both 
suggested by recent observations, raise the possibility of an extended H2 gas phase, well past 
the one inferred by the ^^CO brightness distribution. 

This phase is likely to be warm and diffuse and it may be responsible for the observed high 
total pressures at large galactocentric distances in the Galaxy. In the case of NGC 891 such 
a massive gaseous reservoir may have already been detected through its S(0) line emission. 

We thank the anonymous referee for useful comments that improved the present work. 
P. P. P. would also like to thank Paula Gutierrez Ruiz for many inspiring conversations that 
helped this work through. 
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Fig. 1. — The relative abundance profiles of H2, H, CII, CI and CO for a one-side illuminated 
cloud with total H nuclei density of n = 20 cm~^ (see text). The values of metallicity, FUV 
field and formation rate used are denoted in the upper or bottom left. The effect of an 
enhanced formation rate in making the cloud H2-dominated within much smaller spatial 
scales is apparent. 
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